Slow light in degenerate Fermi gases 
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We investigate the effect of siow light propagating in a degenerate atomic Fermi gas. In particular 
we use slow light with an orbital angular momentum. We present a microscopic theory for the 
interplay between light and matter and show how the slow light can provide an effective magnetic 
field acting on the electrically neutral fermions, a direct analogy of the free electron gas in an uniform 
magnetic field. As an example we illustrate how the corresponding de Haas-van Alphen effect can 
be seen in a gas of neutral atomic fermions. 
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The recent advances in trapping and cooling atoms has 
provided an excellent starting point for studying many 
different types of physical phenomena, ranging from fun- 
damental atomic physics to cosmological aspects [1]. In 
this respect atomic Bose-Einstein condensates have at- 
tracted a lot of interest [2] . Recently several experimen- 
tal groups have succeeded in trapping and cooling also 
fermions [3, 4] well below the Fermi temperature. Fermi 
systems are well known from the study of electron prop- 
erties in materials. Trapped atomic fermions are electri- 
cally neutral and a direct analogy between the magnetic 
properties of these systems and solid state phenomena is 
not necessarily straightforward. We suggest this prob- 
lem can be circumvented if the properties of slow light 
is used, i.e. light with a group velocity as low as meters 
per second [5-7]. The coupling between the slow light 
and the atoms can give rise to some remarkable effects 
such as dragging of the light [8-10] and complete coher- 
ent freezing of the pulse [11-13]. In a similar manner the 
slow light should affect the atomic motion. 

In this Letter we investigate the influence of slow light 
on the mechanical properties of a degenerate Fermi gas 
of atoms. The theory is fully microscopic and based on 
the explicit analysis of the quantum dynamics of atomic 
fermions coupled to the electromagnetic field. In partic- 
ular we use slow light with an orbital angular momentum 
[14, 15]. This allows us to introduce an effective magnetic 
field which acts on the electrically neutral fermions. As 
such we have a typical, often regarded as an academic, 
text book scenario with free electrons moving in a con- 
stant magnetic field. This opens up the possibility to 
study phenomena well known from solid state and con- 
densed matter physics, with all the benefits given by the 
trapped atoms where a range of experimental parameters 
such as atom-atom interactions, particle numbers, the 
shape of the trapping potential etc. can easily be manip- 
ulated. In addition, using light as the effective magnetic 
field is going to be favourable since it is rather difficult 
to control real magnetic fields. As an example we show 
how the de Haas- van Alphen effect is obtained in a neu- 



FIG. 1: a) The level scheme for the Electromagnetically 
Induced Transparency with the probe beam fi p and control 
beam Q c . b) Schematic representation of the experimental 
setup with the two light beams incident on the cloud of atoms. 
The probe beam propagates in the ^-direction. The control 
beam can propagateparallel [12, 13], perpendicular [5] or an- 
tiparallel to the probe beam. 



tral cloud of fermions. Finally we conclude by discussing 
some other aspects and possibilities with slow light in 
degenerate quantum gases. 

Light can be slowed down [5-7] by using the proper- 
ties of Electromagnetically Induced Transparency (EIT) 
[16-19], in which the group velocity of the light, v g , is 
reduced by applying another beam, called the control 
beam. The beams act on the A-type atoms characterized 
by two hyper- fine ground levels 1, and 2, as well as an 
electronic excited level 3, as depicted in Fig. (1). Initially 
the atoms occupy the lowest level 1. 

The atoms are described in terms of the fermionic field- 
operators ^fj (r, t) representing the second-quantized 
wavefunction for the translational motion of atoms in 
the j-th electronic state, with j = 1,2,3. The operator 
tyj (r, t) annihilates an atom positioned at r characterized 
by the internal state j. In what follows, the spatial and 
temporal variables will be kept implicit in (r, t) = 

The atoms interact with two laser beams: A strong 
control laser drives the transition |2) — ► |3), whereas 
a weaker probe field is associated with the transition 
|1) — ► |3) (see Fig.(l)). The control laser has a fre- 
quency u) c , a wave- vector k c , and a Rabi frequency 
f2 c = fl^ exp (zk c .r), where f^°' ) is a slowly varying am- 
plitude. The probe field, on the other hand, is charac- 
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terized by a central frequency u p — ck pi a wave- vector 
k p = k p z, and a Rabi frequency 

n p = ^°) e w+v) (1) 

where Q p ^ is a slowly varying amplitude. In writing 
Eq.(l) we have allowed the probe photons to have an 
orbital angular momentum Hi along the z axis [14, 15]. 

Introducing the slowly-varying atomic field-operators 
$i = * ie iWlt , $ 3 = <5 3e ^i+"p)* and $ 2 = 
fjg'^i+^-^c)^ an d adopting the rotating wave approx- 
imation, one can write the following equations of motion: 

if&i = -^H 1 + V 1 (v)^ 1 + Mi;^s, (2) 

im 3 =( e3 i-^V 2 )$3 + ^3(r)$3 + ^c$2 + ^ P $i (3) 
fr 2 

ih$2 =(e 2 i-2^V 2 )$ 2 + F 2 (r)$ 2 + ^:$ 3 , (4) 

where m is the atomic mass, Vj(r) is the trapping po- 
tential for an atom in the electronic state j, e 21 = 
h (u)2 — uji + lu c — u p ) and £31 = h (lu 3 — u\ — u p ) are, re- 
spectively, the energies of the detuning from the two- and 
single-photon resonances, huj being the electronic energy 
of the atomic level j. 

It is noteworthy that the dissipation of the excited elec- 
tronic state can be included into equation (3) replacing 
£31 by £31 — ^731 and adding the appropriate noise oper- 
ator. The second hyperfme state |2) has usually a small 
decay rate 721 which can therefore be omitted in the cor- 
responding equation (4). Note also that the equations 
of motion (2)-(4) do not accommodate collisions between 
the ground-state atoms. This is legitimate for the degen- 
erate Fermi gas in which s-wave scattering is forbidden 
and only weak p-wave scattering is present [3, 20-22]. 

Since the probe field is much weaker than the control 
field (| fip I <C \Q C \), depletion of the ground-state atoms 
is small. Furthermore, we assume that the two-photon 
detuning e 2 i is sufficiently small. Neglecting the terms 
with $3, V 2 <f>3 and $3 in Eq. (3), one arrives at the 
adiabatic condition [16-19] relating <I> 2 to <&i as: 

* 2 (r,t) = -C*i(r,t). (5) 

where £ = £(r) = Q p /Q c . The condition (5) implies that 
the control and probe beams have driven the atoms to 
the dark state 1 1) — £ |2) representing a special superposi- 
tion between the two hyperfme ground states [16-19]. If 
the atoms are in the dark state, the resonant control and 
probe beams can not populate the upper atomic level 3, 
as the two beams contribute destructively to the absorp- 
tion process. This justifies neglecting the decay of the 
upper atomic level 3 in the equation of motion (3). 

Equation (5) shows that the orbital angular momen- 
tum M of the probe field Q p ~ e lt ^ is transferred into 
the orbital angular momentum of the centre of mass mo- 
tion for atoms occupying level 2. This goes along with 



a general rule saying that the exchange of the orbital 
angular momentum in the electric dipole approximation 
occurs exclusively between the light and the atomic cen- 
tre of mass motion [23]. The rule has been implicitly 
assumed in the initial equations of motion (2)- (4) con- 
taining no contributions due to exchange in the orbital 
angular momentum between the internal atomic states 
and the centre of mass motion. 

Consider now the influence of the slow light on the 
dynamics of the ground state atoms. Using Eqs.(4) and 
(5), one has: 

^(M) = -^(^V 2 + 4-e2 1 - W )( W . 

(6) 

The relationships (2) and (6) provide the following equa- 
tion of motion for the field operator $i(r), 

= ^ [toV + A e// ] 2 $x + y e //(r)$i, (7) 

where 

A e// (r) = ihCV( = -h\(\ 2 VS + z^V|C| 2 (8) 

and 

V eff (v) = V, (r) + (iCf 2 - 2) + hu 21 Id 2 (9) 

are the effective vector and trapping potentials, 

with hu2i = £21 + U 2 (r) — Vi(r). Here the dimension- 
less function ( = e tS fl p ^ / 'Q,^ is characterized by a 
phase S = (k p - k c ) • r + (.(f). In writing Eqs. (7)-(9), 
we made use of the assumprion |£| 2 <C 1. Note that 
such an assumption is not essential in deriving Eqs. (7)- 
(9). By relaxing the condition |(| 2 <C 1, one arrives at 

A-eff = ^ + K| 2 ) C*^C- In such a situation, V e ff 
also experiences modifications. This effect will be ex- 
plored elsewhere. 

It is instructive to note that A e // is generally non- 
Hcrmitian. The Hermitian contribution is due to the 
changes in the phase of £, the non-Hermitian one be- 
ing due to the changes in the amplitude. The non- 
Hermitian part of A e ff can be eliminated by a pseudo 
gauge transformation $1 = $^°- ) exp[— |C| 2 /2], where the 
transformed operator undergoes a unitary evolution. 
Since \(\ 2 <C 1, one can neglect the small changes in the 
amplitude of $1 making the operator A. e ff Hermitian. 
Note also that the probe field Cl p is considered to be an 
incident quantity not affected by the induced motion of 
the ground-state fermions. Consequently the probe field 
fl p undergoes a usual propagation at a group velocity 
v g ~ \Q C \ [16-19] in the ^-direction. 

In this way, we can create an effective vector potential 
through the phase S of the incoming probe beam. The 
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experimental situation is schematically described in Fig. 
(1) where the incoming probe beam is of the form e ^ . 
Suppose that the intensity of the control field does not 
vary considerably within the atomic cloud. If we consider 
co-propagating control and probe beams and choose the 
intensity of the probe beam of the form |O p | 2 ~ r 2 in the 
transversal plane, we obtain the effective vector potential 



is rather intractable. We can, however, obtain analytical 
expressions for the eigenvalues in the limit \£\cto 3> 1, 
where the energies are of the form 



2\£\a (2r 



1) 



(14) 



-M\(\ 2 V(j)-- 



(ye x - xe y ) (10) 



where ao = \(\ 2 R 2 /r 2 is a ratio (typically a < 0.1) be- 
tween the intensities of the probe and the control beam 
at radius R of a cylinder in which the gas is contained. 
Such an external trap can be created by for instance high 
order Bessel beams [24, 25]. It is interesting to note here 
that with this choice of light the effective vector poten- 
tial (10) corresponds to a constant magnetic field in the 
direction opposite to the z-axis, since we have the relation 



with n = 0, 1, 2, ... and q = ... — 2, —1, 0, 1, 2.... This is in- 
deed the Landau result, since for \£\ao » 1 the boundary 
conditions play a less important role. 

With the discrete energy levels given by Eq.(14) we 
are now in a position to calculate thermodynamic poten- 
tials, such as the free energy, as a function of the effective 
magnetic field. The free energy is defined as 



F = Nrj - 



Z(E) 



e {E-n)/kT + i 



dE 



(15) 



B 



eff 



= Vx A eff = - 



2Ma . 



-e z . 



(11) 



The strength of the effective magnetic field is given by 
the orbital angular momentum of the light and can be 
controlled by applying suitable phase and intensity holo- 
grams [15]. It is relatively straightforward to create and 
control high angular momenta of the order of several hun- 
dred £ which consequently controls the effective magnetic 
field. 

Equation (7) describes the effective dynamics of 
trapped noninteracting atoms obeying the Fermi-Dirac 
statistics. Substituting \&i — ^e lEt ' h leads to the follow- 
ing eigenequation for the one particle solution 5', 



where Z(E) is the number of states with energies less 
than E, r/ is the chemical potential, N the number of 
atoms, k the Boltzmann constant and T the temperature. 
In order to proceed we need an expression for the function 
Z(E) using the discrete energy levels in Eq. (14). In 
addition we assume a continuous atomic spectrum in z- 
direction. The resulting function is 



Z{E) 



h 2tt 



(16) 



n,q 



where the sum is over n and q such that Z(E) is real 
and L z is the length of the cloud in z-direction. In the 
limit of £ = S> 1 and with the rescaled temperature 
^ = M we ca l cu l a t e i after carefully examining the sum 
in Eq. (16) [27], the free energy which becomes 



2m 1 



-V 2 + (^)V + 2z( 



too 
R 2 



= E* (12) F = Nr] _ A 



24 



3072 ^ + 35 £ ° 



96 



ete 2 - 



where the external trap V\ (r) is chosen such that V e j / = 
0. After rescaling the radial coordinate, r = xR, and 
using the ansatz * = S,(r)e lq ' t 'e lk '' z we obtain the solution 
in the form of a confluent hypergeometric function 



3 ~ yicos[7re s - \ ] - ^ cos[we i 
(-l) s s 2 sinh[27r 2 6»s] 



(17) 



where e 



+ + d3) 

(E - The result presented in 

Eq.(13) differs from the standard Landau problem [26] 
for free electrons in a homogenous magnetic field in the 
sense that in our case we have to take into account the 
boundary conditions at r = R and the fact that with a 
particular choice of light beam we choose the form of the 
vector potential whereas in the Landau case the vector 
potential is not uniquely defined. It is interesting to note 
at this point the asymmetry induced by the vector poten- 
tial. This is clearly seen in Eqs. (12) and (13) where the 
energy eigenvalue is shifted by 2£aoq- As such, Eq. (13) 



where A = ^^(nn)i, A = and m = 

h 2 \£\ao/2mR 2 . Using the free energy it is straightfor- 
ward to calculate thermodynamic properties such as the 
specific heat 



C =-T 



8 2 F 
9T 2 ' 



(18) 



see Fig. (2), which consequently has an oscillating term 
as well. From Eq.(17) we see that for 9 <C 1 the specific 
heat is linear in 9 whereas the oscillating behavior stems 
from the 1 / ft£l dependence in the cos-term. 

In order to have a dominating oscillating part 9 must 
not be too large, < 1, otherwise the oscillating term is 
exponentially damped. With present experimental cool- 
ing and trapping techniques temperatures of the order 
of 9/eo = T/Tp <~ 0.1 are readily achievable. Hence an 
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FIG. 2: a) The specific heat as a function of the rescaled 
temperature 6. The lowest line corresponds to eo = 1, middle 
one to £o = 10 and the top one to £o = 30. b) The specific heat 
as a function of the inverse effective magnetic field, eo = r//hQ 
where the three curves correspond to 6 = 0.1, 0.5 and 1.0 (top 
one). The inset shows a magnification of the region for small 
eo and for 6 = 0.1. 



eo = rj/Ml of the order of one would be preferable. For 
a homogeneous cloud in a cylindrical trap e is given by 
£o = {Nj^Stv) 3 /(ao£) which can become small if one 
considers a large aspect ratio trap and taking into ac- 
count the fact that the term aal can reach large values 
of the order of 100. In Fig. (2) we show the specific 
heat calculated from Eq. (18) where C is in units of 
fc VgmM j n ca i cu i a ti n g C we used the exact equation 
(15) for F to show the full range of the £o-dcpcndcncc. 

In this Letter we have shown how light with an orbital 
angular momentum can be used to create an effective 
magnetic field in a degenerate gas of electrically neutral 
atomic fcrmions. As an example on how the slow light 
can be used we calculated the free energy for the trapped 
degenerate fermions and found that the atomic gas shows 
a de Haas- van Alphen type behavior where oscillations in 
thermodynamic properties depend on the inverse effec- 
tive magnetic field strength. There are other intriguing 
phenomena such as the quantum Hall effect which can be 
studied using cold fermionic gases and slow light with an 
angular momentum. In addition, if the collisional inter- 
action between the atoms is taken into account slow light 
can be used to study the magnetic properties of a super- 
fluid atomic Fermi gas [28]. Recent advances in spatial 
light modulator technology enables us to consider rather 
exotic lightbeams [29]. This will allow us to study the 
effect of different forms of vector potentials in quantum 
gases. In particular the combined dynamical system of 
light and matter could give important insight into gauge 
theories in general. It is certainly tempting to push the 
analogy further and study phenomena from high energy 
physics in ultracold samples of atoms. 
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